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Abstract 

We investigate the spatiotemporal dynamics of coupled circle map 
lattices, evolving under synchronous (parallel) updating on one hand 
Q \ and asynchronous (random) updating rules on the other. Synchronous 

' evolution of extended spatiotemporal systems, such as coupled circle 

j-H . map lattices, commonly yields multiple co-existing attractors, giving 

rise to phenomena strongly dependent on the initial lattice. By marked 
contrast numerical evidence here strongly indicates that asynchronous 
evolution eliminates most of the attractor states arising from special 
, sets of initial conditions in synchronous systems, and tends to yield 

^ I more global attractors. Thus the phenomenology arising from asyn- 

' chronous evolution is more generic and robust in that it is obtained 

, from many different classes of initial states. Further we show that in 

\ parameter regions where both asynchronous and synchronous evolu- 

' tion yield spatio-temporal intermittency, asynchronicity leads to bet- 

. ter scaling behaviour. 
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It is well-known that spatially extended systems undergoing temporal 
evolution show the presence of a large number of multiple co-existing attrac- 
tors In the case of coupled map lattices, which constitute simple models 
of spatially extended systems p|, |[ evolving under synchronous up-dates, it 
has been seen in numerous examples that system attractors shows strong 
sensitivity to initial conditions even in parameter regimes where the evolv- 
ing dynamics is spatio-temporally periodic in nature 0. This multiattractor 
property has significant consequences for problems like control and synchro- 
nisation 1^ in physical, chemical, biological and engineering contexts. The 
stability of such attractors to perturbations and noise have also been studied 
in the case of globally coupled systems 0. 

A typical CML consists of dynamical elements on a lattice which interact 
with suitably chosen sets of other elements, and evolve via synchronous or 
parallel updates wherein the dynamical elements at each lattice site are up- 
dated simultaneously. However, there have been several attempts to study 
CML-s which evolve via asynchronous evolution, that is, one in which the 
updates at lattice sites are not concurrent, but sequential instead. The study 
of asynchronous updates is considered to be interesting for several reasons. 
Notably, neurophysiological systems like neurons and neuron groups evolve 
asynchronously, and lattice dynamical models of such phenomena must of 
course employ asynchronous updating schemes. Therefore it is of importance 
to investigate the effects of asynchronicity in prototype models [6-13]. These 
effects can be quite significant, e.g. it has been argued that asynchronous 
up-dates can alter the universality class of spatio-temporal intermittency [^]. 

In this paper, we observe that asynchronicity can have a very important 
physical effect. Systems which show the existence of multiple co-existing 
attractors which depend on the state of the initial lattice when evolved via 
synchronous updates, evolve to the same global attractor from many kinds of 
initial states under asynchronous evolution. Thus asynchronicity wipes out 
multiple co-existing attractors and leads to a generic global attractor which 
is robust to the evolution of different classes of initial conditions. 

We demonstrate this result for the following specific example. The space 
on which our CML dynamics occurs is a discrete 1-dimensional chain, with 
periodic boundary conditions. The sites are denoted by integers i, i = 
1, . . . ,N, where N is the linear size of the lattice. On each site is defined a 
continuous state variable denoted by Xn{i), which corresponds to the physical 
variable of interest, with n denoting the discrete time. Here the local on-site 
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map is chosen to be the sine circle map, a map which has generated much 
research interest: 

K 

f{x) = X -\- Lo sin(27ra;) 

27r 

where < x < 1 . The parameter K indicates the strength of the nonhnearity, 
and is chosen to be 1 here. 

These local maps are coupled through their nearest neighbours, and the 
coupling form is the discretized Laplacian form (i.e. future coupled) [Q], 
with the strength of coupling given by e. Now extensive results have been 
obtained for the standard parallel time evolution of lattices of circle maps, 
i.e. the scheme in which all individual maps of the lattice are iterated forward 
simultaneously Q]. This synchronous evolution implies the dynamics: 

= (1 - e)/(x„W) + |{/(x„(z - 1)) + f{xn{i + 1))} mod 1 (1) 

It is known 0, || that varying e and cu, under the usual parallel undat- 
ing, yields various dynamical phases, such as: the synchronized fixed point, 
i.e. spatial period 1 temporal period 1 (SlTl), spatio-temporally periodic 
solutions e.g. spatial period 2 temporal period 1 (S2T1) and spatial period 
2 temporal period 2 (S2T2), spatial intermittency, spatiotemporal intermit- 
tency etc. For parallel updating an important feature of the spatiotemporal 
dynamics is the existence of multiple co-existing attractors. This implies 
that the resultant attractor is very strongly dependent on the spatiotempo- 
ral features of the initial lattice. Thus synchronous evolution can yield many 
different behaviours under different initial preparations even for identical pa- 
rameter values. Detailed phase diagrams of the behaviour of this CML under 



synchronous evolution have been obtained [Q, . 

Here we focus on randomly updated CMLs, i.e. a CML where the el- 
ements of the lattice do not update simultaneously, but update one after 
another in random sequence. Time now is measured in units of one complete 
updating sweep of the lattice (in random order), i.e. after one unit of time, all 
sites in the lattice have been updated. This random updating scheme should 
help us test the implications of asynchronicity in extended systems, and help 
us gauge the degree of robustness of various physical features emerging under 
conventional synchronous evolution. 

We study the different e and u regions of this coupled circle map system, 
under three classes of initial lattices: (i) spatial period 2 initial condition; (ii) 
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spatial period 2 initial lattice, with a kink; and (iii) random initial lattice. We 
observe the spatiotemporal dynamics under parallel updating on one hand 
and completely asynchronous updating on the other, and find out how this 
"phase diagram" changes under asynchronicity, for the different classes of 
initial states mentioned above. 



Case I: In the region with very high coupling parameter e and low for 
instance around e ~ 0.97, u ~ 0.02, parallel updating yields spatial period 
2, temporal period 2 behaviour when evolved from spatial period 2 initial 
lattices. But for the case of the initial periodic lattice having a kink, and for 
the case of random initial lattices, spatio-temporal intermittency is obtained. 
(See for a detailed phase diagram). 

By contrast, under asynchronous updating in this parameter regime one 
obtains spatiotemporal fixed points for all initial conditions: spatial period 
2, spatial period 2 with a kink, and random initial lattices. This fixed point 
is at ^ = f{6) ~ u for low u. (See Fig. 1.) Note that this is in agreement 
with earlier observations that under strong coupling asynchronicity has the 
effect of regularising the system |]12| . 

Thus, under randomness in updating rules, spatiotemporal dynamics is 
no longer sensitive to initial conditions. The spatiotemporal fixed point is 
now the global attractor of the system. 



Case II: In the region around e ~ 0.8, u ~ 0.2, asynchronicity yields spa- 
tiotemporal chaos from all initial conditions. By contrast, parallel updating 
yields spatiotemporal periodicity when evolved from spatial period 2 initial 
lattices (with or without kink), and yields spatiotemporal intermittency when 
evolved from random initial lattices (see Fig. 2). 

Also note here, that while asynchronicity had the effect of regularising 
the system under strong coupling (as seen above in Case I), under weaker 
coupling asynchronicity has the effect of inducing spatiotemporal disorder. 
In this example for instance, we see that parallel updating gives rise to spa- 
tiotemporal periodicity for certain classes of initial states, while asynchronous 
evolution always leads to spatiotemporal chaos. 



Case III: Around the region e ~ 0.8, u ~ 0.08, asynchronicity yields 
spatiotemporal intermittency for all initial conditions. In contrast, parallel 
updates yield exact spatiotemporal periodicity for period 2 initial lattices. 
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and spatiotemporal intermittency for random lattices. 

Interestingly, when both asynchronous and synchronous evolution yield 
spatiotemporal intermittency from random initial lattices, the scaling expo- 
nents obtained from the probability distribution of laminar lengths in space 
and time are quite distinct for the two cases. The intermittent dynamics 
arising from asynchronous evolution exhibits scaling in both time and space, 
with the spatial and temporal scaling exponents being approximately the 
same (0 ~ 3). Synchronous evolution on the other hand leads only to tem- 
poral scahng, and not good spatial scahng (see Fig. 3). 

We note that similar phenomena arc observed for coupled logistic map lat- 
tices as well, i.e. a CML with the local on-site map being f{x) = rx(l — x), 
with r — A. See figures 6 and 7 for two examples, one for high coupling 
strength e = 0.9 and the other for low coupling strength e = 0.1. Clearly 
here too asynchronous evolution is insensitive to differences in initial con- 
ditions, whereas multiple attractors co-exist for the case of the traditional 
synchronous evolution. 

In summary then, we have investigated the spatiotemporal dynamics of 
coupled map lattices evolving under asynchronous updating rules. We have 
shown that asynchronicity has a pronounced effect on the spatiotemporal dy- 
namics: it yields a more global attractor where multiple attractors co-existed 
for the case of parallel updating. In that sense asynchronous updating yields 
more generic and robust phenomena in extended systems. Conversely, intro- 
ducing some degree of asynchronicity in an extended system may help lead 
the system to its most generic attractor from any initial condition, however 
special. The nature of the generic attractor i.e. whether regular or disor- 
dered, appears to depend on the strength of the coupling, as well as the 
synchronicity /asynchronicity. In the case of spatio-temporal intermittency, 
asynchronicity leads to better spatial scahng. Thus, synchronous and asyn- 
chronous updates can lead to different universality classes of spatio-temporal 
behaviour for systems which are otherwise identical in all respects. There- 
fore asynchronicity constitutes a relevant perturbation in the evolution of 
extended systems, and can perhaps be used to direct spatially extended sys- 
tems to desired global attractors. We hope our observations will be useful 
towards the understanding the diversity of phenomena which can be seen in 
spatially extended systems and for controlling their behaviour. 
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Figure 1: Area of e-u parameter space occupied by the spatiotemporal fixed 
point for the case of parallel updating (top panel) and completely asyn- 
chronous updating (bottom panel), for different initial conditions: (a) period 
2 initial lattice (b) period 2 with a kink and (c) random initial conditions. 
Here lattice size = 100. 



Figure 2: Space-time density plots showing the evolution of a circle map 
lattice of size = 100 over 100 iterations via parallel updating and via 
completely asynchronous updating. Each slice parallel to the x axis displays a 
snap shot of the spatial profile at some instant of time. Here e = 0.8, u = 0.2. 
The lattices are evolved from (top to bottom) period 2 initial lattice, period 
2 lattice with a kink at i = 37 and random initial lattice. 
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Figure 3: Probability P of occurence of (a) temporal laminar regions and 
(b) spatial laminar regions vs length of laminar regions / (on a log- log plot), 
for the case of asynchronous updating (solid squares) and parallel updating 
(solid triangles). The solid hnes are the best fit straight hnes, with slope 
equal to —3.0 for asynchronous updating and —2.4 for parallel updating in 
(a), and with slope —3.2 for asynchronous updating in (b). These results are 
for circle map lattices of size = 1000, with e = 0.8, a; = 0.08, evolved from 
random initial conditions. Sample size is ~ lO'^. 
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Figure 4: Space-time density plots showing the evolution of a logistic map 
lattice of size N = 100 over 100 iterations via parallel updating and via 
completely asynchronous updating. Each slice parallel to the x axis displays 
a snap shot of the spatial profile at some instant of time. Here e = 0.9. The 
lattices are evolved from (top to bottom) uniform (spatial period 1) initial 
lattice, spatial period 2 initial lattice and random initial lattice. 



Figure 5: Space-time density plots showing the evolution of a logistic map 
lattice of size N = 100 over 100 iterations via parallel updating and via 
completely asynchronous updating. Each slice parallel to the x axis displays 
a snap shot of the spatial profile at some instant of time. Here e = 0.1. The 
lattices are evolved from (top to bottom) uniform (spatial period 1) initial 
lattice, spatial period 2 initial lattice and random initial lattice. 



